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ABSTRACT 

We report new results from recent Combined Array for Research in Millimeter-wave Astronomy 
(CARMA) observations of both continuum and HCO+(1-0) line emission at A 3.4 mm from W3-SE, 
a molecular core of intermediate mass, together with the observations of continuum emission at A 1.1 
and 0.85/0.45 mm with the Submillimeter Array (SMA) and the James Clerk Maxwell Telescope 
(JCMT), respectively. A continuum emission core elongated from SE to NW, with a size of ~10", has 
been observed at the millimeter and submillimeter wavelengths. The dust core has been resolved into 
a double source with the SMA at A 1.1 mm. The angular separation between the two components 
is ~4". Together with the measurements from the Spitzer Space Telescope and the Midcourse Space 
Experiment (MSX) at mid-IR wavelengths, we determined the spectral energy distribution (SED) of 
the continuum emission from W3-SE and fit it with a thermal dust emission model. Our best fitting of 
the SED suggests the presence of two dust components with different temperatures. The emission at 
mm/submm wavelengths is dominated by a major component which is characterized by a temperature 
of Td=:41±6 K with a mass of 65±10 Mq. In addition, there is a weaker hot component 400 K) 
which accounts for emission in the mid-IR, suggesting that a small fraction of dust has been heated 
by newly formed stars. We also imaged the molecular core in the HCO"'"(1-0) line using CARMA at 
an angular resolution '^6". In the central region of ~50", the integrated IICO"'"(1-0) line emission 
shows a main component A that coincides with the dust core, as well as two sub-structures B and 
C which are located N and SE of the dust core, respectively. With the CARMA observations, we 
have verified the presence of a blue-dominated double peak profile toward this core. The line profile 
cannot be explained by infall alone. The broad velocity wings of the line profile suggest that other 
kinematics such as outflows within the central 6" of the core likely dominate the resulting spectrum. 
The kinematics of the sub-structures of B and C suggest that the molecular gas outside the main 
component A appears to be dominated by the bipolar outflow originated from the dust core with a 
dynamical age of >3xl0^ yr. Our analysis, based on the observations at wavelengths from millimeter, 
submillimeter to mid-IR suggests that the molecular core W3-SE hosts a group of newly formed young 
stars and protostars. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - ISM: molecules - ISM: jets and 
outflows - radio lines: ISM - stars: formation 



1. INTRODUCTION 

Among the dynamical processes in star formation re- 
gions, infall appears to be critical in understanding the 
processes of mass accumulation in star formation cores 
and the growth of protostars. The protostellar infall 
rate depends on the surface density of the protostellar 
core, which usually can be determined from observations 
of self-absorptions of the embedded molecular clumps in 
it. However, it is difhcult to unambiguously identify in- 
fall based on low angular-resolution observations because 
the spectral features of a complex with bipolar molecu- 
lar outflows, rotating structures or multiple cores can 
produce similar characteristics to those expected from 
self-absorption in infalling clumps. In many young star 
formation regions that show spectral signatures of in- 
fall from single-dish observations, the single-dish spec- 
tral features often turn out to be dominated by outflows 
and other kinematic features with further high-resolution 
observations using interferometer arrays. 
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W3-SE is a den se molecular core located 2 kpc away 
(|Blitz et al.lll982D . SE of the well-known massive star 
formation region W3-Main . The W3-SE core was first 
imaged by iTieftrunk et all (|1998f ) in the NH3{1,1) and 
(2,2) lines. The low kin etic temperature of 25 K based 
on the NH3 line ratio (ITieft ru nk et al.l Il998f) and the 
absence of IRAS emission ( Jiiina et al.l [1999') implied 
that this core is probably in an early evolutionary stage. 
A virial mass of ^^300 Mq was inferred from the line 
widths, indicating that this core may be a birthplace 
for intermediate- or high-mass stars. However, to date 
there has been no evidence for OB stars present in the 
core because of the non-detection from th e radio contin- 
uum emission at 6 cm with the VLA (Ticft runk et al.l 
119981 ). On the other hand, a possible infall signature 
in W3-SE was suggested by IWu et all (|2007D based on 
observations of multiple molecular lines with the IRAM 
30-m telescope. Both the HCO+ (1-0) and (3-2) hues 
show an asymmetric double profile characterized by a 
primary peak blue-shifted with respect to a peak veloc- 
ity of the optically thin lines such as C^^O(l-O), while a 
secondary peak is red-shifted. The possible signature of 
infall in W3-SE suggested from the single-dish observa- 
tions needs further verificati on with impr oved, high reso- 
lution observations. Besides. [Ruch et aD |2007 ) reported 
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a detection of a mid-IR point source in the region of W3- 
SE and determined the flux densities in the four bands of 
the Infrared Array Camera (IRAC) on the Spitzer Space 
Telescope. The color-color diagram in their paper indi- 
cates that this source is most likely to be a YSO (Class 
I or earlier) with large excess emissions at longer wave- 
lengths. 

In this paper, we present new results based on ob- 
servations at A 3.4, 1.1, 0.85/0.45 mm with the Com- 
bined Array for Research in Millimeter-wave Astronomy 
(CARMA), the SubmiUimeter Array (SMAj!| and the 
James Clerk Maxwell Telescope (JCMTj^, respectively, 
along with the data at mid-IR wavelengths obtained from 
both the Spitzer Space Telescope and Midcourse Space 
Experiment (MSX) archives. §2 describes observations 
and data reductions, §3 presents results from the anal- 
ysis of the continuum emission at multiple wavelengths 
and the CARMA observations of the IICO+(1-0) hue, 
and §4 discusses the results and their implications for 
astrophysics. Finally, §5 gives the summary and conclu- 
sions. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. CARMA Observations at X 34 mm 

W3-SE was observed at A 3.4 mm with the CARMA 
telescope in the D-configuration on June 24th and June 
27th, 2008, with an angular resolution of ~ 6". The ob- 
servations were centered at RA(J2000) = 02'^25'"54^50, 
DEC(J2000) = 62°04'11".0, which is used as the refer- 
ence position in the following analysis. Two broad bands 
with BW of 500 MHz were set at 86 and 89 GHz to image 
the continuum emission. Two USB narrow bands with 
BW of 8 and 32 MHz were tuned to the HCO+(1-0) line 
at 89.188 GHz, and the two LSB narrow bands were used 
for the SiO(2-l) line at 86.243 GHz. The highest veloc- 
ity resolution is 0.42 km s~^ for the 8-MHz narrow-band 
data with 63 spectral channels. The observations were 
Doppler-tracked at ^lsr = -38.9 km s~^. 

The observations include two tracks. The first one 
lasted for 5 hours using 14 antennas with a 19-point 
mosaic, and the second one was taken for 4 hours us- 
ing 13 antennas with a 7-point mosaic. Because the first 
track was on a maintenance day, several antennas were 
not used in parts of the track. The corrupted uv data 
were flagged during the off-line data reduction. The sec- 
ond track was carried out in good weather, and the data 
quality is excellent with Tgys between 150 and 200 K. 

All the data reduction s and imaging pr ocesses were 
carried out in MIRIAD (|Sault et al.lll995D . The data 
from the two tracks were calibrated separately with the 
calibrators listed in Table 1. The bandpass and complex 
gain calibrations were processed following the standard 
procedure for the CARMA data reduction. Flux density 
scales were calibrated using Uranus (7.2 Jy) and MWC 
349 (1.3 Jy) at A 3.4 mm for the two tracks, respectively. 

^ The SubmiUimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 

The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Particle Physics and Astronomy 
Research Council of the United Kingdom, the Netherlands Organ- 
isation for Scientific Research, and the National Research Council 
of Canada. 



After calibration the two datasets were merged and im- 
ages were made with a linear mosaic algorithm using the 
program INVERT. 

2.2. Combining CARMA Data with IRAM 30-m Data 

The shortest baseline of the CARMA observations is 
2.7 kA. To eliminate the effects of missing short spac- 
ing, we combined the CARMA and IRAM 30- m telescope 
data which cover the central 45"x45" region (jWu et aTl 
120071 ). We used the MIRIAD program IMMERGE to 
linearly combine the single-dish image and the cleaned 
CARMA image. The uv ranges sampled by the two ob- 
servations are well overlapped, which are used to deter- 
mine the scaling factor for the flux density of the single- 
dish data. The amplitude of the IRAM 30-m data was 
re-scaled by applying the scaling factor determined by 
comparing the amplitude values of the overlapped visi- 
bilities in the FFT domain. There is a boundary effect 
on the edge of the combined image of 45"x45" due to 
limited spatial coverage of the IRAM observations. The 
effect is negligible in our analysis of the line emission in 
the central 20"x20" region. 

2.3. SMA Observations at X 1.1 mm 

We also made observations at A 1.1 mm with the SMA 
in the compact array with an angular resolution of ^^2.5". 
The SMA observations were centered at the same posi- 
tion as that used in the CARMA observations. The fre- 
quencies were tuned to 266.7/276.7 GHz for LSB/USB 
with bandwidths of 2 GHz for each of the side bands. 
The typical Tgys in the observations is 180-200 K. 

The data reduction and imaging were made using 
MIRIAD. The bandpass shapes and complex gains were 
calibrated using QSOs 3C454.3 and 0244+624, respec- 
tively. The flux density scale was bootstrapped from 
Uranus using the SMA planet model. The continuum 
image of W3-SE was constructed by combining the visi- 
bility data from the two side bands with a rms noise of 
1.5 mJy. 

2.4. JCMT Observations at X 0.85/0.45 mm 

The W3-SE region was observed using JCMT with 
SCUBA at A 0.85 and 0.45 mm on September 15, 2004 
(PID M04BC17). The data were calibrated and com- 
bined into intensity images with the ORAC-DR pipeline. 
The FWHM beam sizes of the JCMT images are 14.5" 
and 7.5" at A 0.85 and 0.45 mm, respectively. The JCMT 
images were re-grided with the program REGRID in 
MIRIAD to align the coordinate axes with those in the 
CARMA images. 

2.5. Spitzer and MSX Data at Mid-IR 

The entire region of W3 was observed by R. Gchrz 
(the PI in GTO program PID 127) with the IRAC 
iFazio et al.ll2004D on January 10, 2004 an d the Multi- 
band Imaging Photometer for Spitzer fMIPS: [Rieke et aTl 
I2OOI on January 10, 2004 on board the S pitzer Space 
telesco pe. The results have been published bv lRuch et al.l 
(|2007D . and a point source which is detected in all four 
IRAC bands in the W3-SE region was reported. In order 
to search for additional weaker mid-IR sources to identify 
the possible counterparts of the radio sources that have 
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been detected at mm/submni wavelengths, we have fur- 
ther investigated the W3-SE region by re-reducing the 
archive data of IRAC and MIPS (GTO program PID 
127). 

The basic cahbrated data (BCD) products of IRAC 
are provided by the Spitzer Science Center (SCC) with 
the IRAC pipeline version S18.7.0. A post-BCD script 
hdr_mask.pl provided by SSC was used to fix the poten- 
tial saturation problem. The background correction and 
mosaic were made with the version 1 8.2.0 of the post- 
BCD reduction software MOPEX (IMakovoz fc KhanI 
|2603 ). producing a mosaic image with 1.22" per pixel. 
The point-source extraction and the aperture photome- 
try were done with the package APEX in MOPEX. Point- 
spread function (PSF) fittings were carried out in the 
process of detecting point-like sources, and those with 
signal-to-noise ratio higher than 6 are identified as point 
sources. An aperture with radius of 2.4" and an annu- 
lus between 2.4" and 7.3" were used for the aperture 
photometry to determine the flux densities of the point 
sources. Aperture correction, color correction and pixel- 
phase correction were carried out according to the pro- 
cedures described in the IRAC Data Handbook (version 
3.0). In addition, the flux density of the W3-SE core was 
determined using the procedures of aperture photometry 
for extended sources with a radius of 8.5" in order to com- 
pare with the flux densities of millimeter/submm in the 
same region and to determine the SED. We noted that 
it was difficult to make a reliable annulus to represent 
the background levels due to the background/foreground 
stars. Therefore, to avoid obvious stellar emission, we 
measured fiux densities in several regions which are close 
to W3-SE to calculate the background levels. An overall 
fractional error of ~15% was assessed for the photome- 
try, dominated by the uncertainties in the corrections for 
the background/foreground emission levels (~10%) and 
the aperture correction (^10%). 

The BCD products of MIPS were pre-processed by the 
MIPS pipehne version S16.1.0. The W3-SE region is 
badly saturated in the 70 and 160 /xm images due to the 
strong emission of W3-Main; thus only the 24-/im data 
were used. A post-BCD script flatfield.pl provided by 
SSC was used to make an extra flat-fielding to remove 
the potential residual background gradient, and then 
the background correction and mosaic were made with 
MOPEX. The aperture photometry was carried out with 
the program aperture in MOPEX with an aperture ra- 
dius of 13.5", and color/ aperture corrections were made 
following the MIPS Data Handbook (version 3.3.1). An 
uncertainty of 15% in the flux densities is dominated by 
the errors in corrections for the background/foreground 
emission (^^10%) and the color/ aperture corrections. 

In addition, we also investigated the MSX images of 
W3-SE in all the four bands (8.3, 16.7, 12.1 and 21.3 
/zm) with an angular resolution of ~20". We re-grided all 
of the Spitzer and MSX images in MIRIAD to align the 
coordinate axes with those in the CARMA images. The 
photometry was carried out with the similar methods 
as those for the IRAC images. The uncertainties of the 
photometry on the MSX images are large ('^40%) due 
to poor sensitivities, low angular resolutions and severe 
contamination from background/foreground emission. 

3. RESULTS 



3.1. Continuum Emission 
3.1.1. CARMA A 3.4 mm 

Figure 1(a) shows the CARMA image of W3-SE at 
A 3.4 mm, revealing an un-resolved dust core plus a tail- 
like structure extended in SE. We determined a peak po- 
sition ((5(R.A.) = -5.5±0.2" and (5(DEC) = -0.9±0.2") 
and peak intensity (24±2 mJy beam~^) for the contin- 
uum emission from the dust core based on a Gaussian 
fitting to the data. Integrated over the core region of 
15"xl5", a total flux density of 42±10 mJy was deter- 
mined. 

3.1.2. SMA \ 1.1 mm 

Figure 1(b) shows the SMA image at A 1.1 mm. The 
dust core has been resolved into two continuum emission 
components (SMA-1 and -2), and the overall morphology 
of the W3-SE dust core shows an elongation from NW to 
SE. The primary (SMA-1) and secondary (SMA-2) com- 
ponents of the double source are located NW and SE of 
the CARMA peak position. From flts with two Gaus- 
sian components and a linear background, we obtained 
peak intensities of 0.21±0.01 and 0.17±0.01 Jy beam"! 
for SMA-1 and SMA-2, respectively. The angular sepa- 
ration of the double components is ~4". The total flux 
densities and sizes of the double components are listed in 
Table. 2. The peak position of the emission at A 3.4 mm 
observed with CARMA agrees well with the centroid po- 
sition of the double components observed with the SMA. 
Wc integrated the continuum flux density from the core 
in the SMA 1.1-mm image using IMAGE in AIPS by 
adding a zero spacing flux density of 2 Jy which was de- 
rived from multiple Gaussian fits to the visibility data, 
and determined that the total continuum flux density 
at 1.1 mm from the W3-SE core (-15"xl5") is 1.4±0.4 
Jy. The large uncertainty in the flux density includes all 
possible error sources, dominated by the missing short 
spacings in the SMA observations. 

3.1.3. ,JCMT A 0.45/0.85 mm 

Figures 1(c) and 1(d) show the continuum images ob- 
served with JCMT at A 0.85 and 0.45 mm, respectively. 
With the angular resolutions 14.5" and 7.5" at A 0.85 and 
0.45 mm, the emission at submm wavelengths from W3- 
SE has been slightly resolved, showing an extended, ~ 1' 
halo, which is better seen in the JCMT image at A 0.85 
mm (Figure 1(c)). An elongated core, 15", is clearly 
shown at A 0.45 mm in Fig. Id. The core-halo morphol- 
ogy of the continuum emission observed at submm with 
the JCMT is consistent with CARMA and SMA obser- 
vations at A 3.4 and 1.1 mm, suggesting that the double 
compact continuum core (SMA-1 and SMA-2) revealed 
by the SMA is encompassed by a halo. Gaussian fits 
to the JCMT continuum images give peak positions (- 
8.2"±0.4", 1.2"±0.6") and (-8.8"±0.4", 0.5"±0.4") for 
the A 0.45 and 0.85 mm emission, respectively. Within 
the JCMT beams, these two positions agree with the po- 
sition of the continuum peak observed with CARMA at 
A 3.4 mm and the centroid position of the double compo- 
nents derived from the SMA observations at A 1.1 mm. 

From the JCMT images, integrating over emission re- 
gion with signals greater than 5 cr, we derived a total 
fiux density 8±1 and 75±8 Jy at A 0.85 and 0.45 mm, 
respectively. All the measurements of fiux densities at 
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millimeter and submm wavelengths are summarized in 
Table 2. 

3.1.4. MIPS and IRAC at Mid-Infrared 

In the MIPS image at 24 /im (Figure 2(a)), the direc- 
tion of the elongation in the 24-/im emission is in NW- 
SE (P.A.«-25°). In comparison with the distribution of 
the emission at millimeters/submm, we noticed that in 
the region located E of the SMA double sources there 
is a large excess emission at 24 /im, indicating a dis- 
placement between the hot and cold dust emission in the 
W3-SE core. 

Based on the IRAC images, a group of discrete emis- 
sion sources were revealed in the core of W3-SE (see the 
3.6 fim image in Figure 2(b)). By fitting PSF, three point 
sources were detected in the central core region of W3 - 
SE. One of them has been reported in lRuch et al.l ()2007l ) 
as the IRAC point source 56 (hereafter IRAC 56), the 
only source that can be fitted as a point source in all 
of the four IRAC bands. The flux densities of IRAC 56 
determin ed from this pape r are in agreement with the 
results of IRuch et all (l2007l ). Located ~4" NW to IRAC 
56, a weak point source (hereafter IRAC 56a) was de- 
tected at 4.5, 5.8 and 8.0 /im, corresponding to a weak 
emission feature at 3.6 /xm which failed to be fitted as 
a point source. The detection of this source is also in 
good agreement with the study of Ruch et al. (Ruch 
2009, private communication). This source appears to 
coincide with the core SMA-1, showing a large excesses 
in the emission at the longer wavelengths (5.8 and 8.0 
/xm) as compared with that at 4.5 and 3.6 /im. The 
large excesses in emission at longer wavelengths suggest 
that probably a protostar is embedded in SMA-1, and 
extinctions at shorter wavelengths become considerable. 
Located ~2" SW to IRAC 56 there appears to be an ad- 
ditional point source (hereafter IRAC 56b) detected at 
4.5 /im, but we failed to fit it as a point source due to the 
contamination from the profile wing of the strong point 
source IRAC 56 and possible extended PAH or H2 emis- 
sion. The nature of this source remains uncertain. In 
addition, we note that there appears no significant emis- 
sion at the mid-IR wavelengths toward the cold dust core 
SMA-2. 

Following the procedures of the aperture photometry 
described early, we determined the mid-IR flux densities 
from the core region of W3-SE. The total flux densities of 
W3-SE at mid-IR wavelengths from the IRAC and MIPS 
images, as well as the values from the MSX images, are 
summarized in Table 3. 

3.1.5. Dust Content and SED 

Figure 3 shows the spectral energy distribution (SED) 
of W3-SE covering the wavelengths from millimeter to 
mid-IR. We fit the SED with a thermal dust model to 
determine the physical parameters of the dust compo- 
nents in the core of W3-SE. Given the dust temperature 
Td and the total mass of gas and dust Mtotai, the con- 
tinuum emission at the frequency v from thermal dust 
in a molecul ar core subtending a solid angle Vl can be 
described as ()Wilson et al.ll2009[ ) 

S, = B,{n){\ - e-^-)n (1) 

where 



In the optically thin case < 1, Equation (1) can be 
approximately expressed as 

where Ki, is the dust opacity per unit dust mass, g is the 
density ratio of II2 gas to dust, D is the distance, and 
Bu{Td) is the Planck function with a dust temperature of 
Td. A dopting k(300 GHz) = 1.4 cm ^ g~^ as the reference 
value (jOssenkopf fc HenninglllQQ^ . we assumed a power 
law of = kq(v/vqY ()Hildebrandl I19831) . where /3 is 
the power-law index for dust opacity. In addition, we 
adopted g = 100 and D = 2 kpc in the following analysis. 
Thus, for each dust component, the thermal dust model 
has three free parameters (Mtotai, Td and /?). 

Figure 3(a) shows least square (LSQ) fitting to the 
SED with a total of thirteen measurements from submm 
to mid-IR on the assumption of optically thin dust emis- 
sion. We found that a model with two thermal dust 
components appears to adequately fit the data. The six 
free physical parameters of the model can be determined 
by the LSQ fitting to the data. Our best fitting suggests 
that a cold, massive dust component accounts for the 
continuum emission at the millimeters/submm as well as 
a large fraction of the continuum emission at 24 fim. A 
total mass of Mtotai = 70±7 Mq, a dust temperature of 
Td — 35±1 K, and a dust opacity index of /3 = 2.0 ±0.1 
were derived from the fitting for this cold, massive dust 
component. In addition, the emission bump at mid-IR 
requires a secondary dust component with a small mass 
but a high temperature. Our fitting gives Td = 4OOI5Q 

K, 13 = O.btoi and A'/totai < 5 x lO"'^ Mq for the sec- 
ondary component. We note that the uncertainties of 
the derived parameters for the secondary component are 
large due to the large uncertainties in the mid-IR mea- 
surements, especially the ones from MSX. 

Although the optically thin approach is valid for the 
emission from the cold dust component at millimeter and 
submm wavelengths and the emission from the hot com- 
ponent at mid-IR, the optical depth of the cold compo- 
nent may become significantly large at the wavelengths 
of mid-IR. On the assumption that the steep power law 
of the dust opacity {f3 — 2) holds up to the mid-IR, the 
mean optical depth of r^y ^160 at A 24 /im for the cold 
dust in the region of the double source is inferred using 
the values of Td, and Aftotai derived from the optically 
thin approach. Thus, the optical-depth effect can not be 
ignored for the dust emission at the shorter wavelengths. 
Figure 3(b) shows a good fit to the observed SED data, 
considering self-absorption in the cold dust with Equa- 
tion (1) and assuming that cold component docs not over- 
lap the hot dust emission region. Consequently, in the 
self-absorption model the turnover frequency i't=i, where 
r(z^) — 1, is shifted to the higher frequency of 1070 GHz 
as compared with the value (j^t=i = 990 GHz) derived in 
the optically thin approach. The self-absorption model 
gives a 15% decrease in Mtotai and a 25% increase in Td, 
i.e., Td = 47 K, and Mtotai = 60 M©, if the power-law 
index (3 — 2 remains unchanged at mid-IR. In the case 
that the power-law of the dust opacity becomes flat, e.g., 
(3 = 0.5 at the frequencies above 1000 GHz, then 
at A 24 /im and the optical depth effect at mid-IR be- 
comes not so critical. In any case, the uncertainties in 
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the determination of Td and Mtotai appears to be mainly 
due to the uncertainty in the dependence of the opacity 
on frequency at mid-IR. Thus, we take the mean values 
of the physical parameters derived from the two extreme 
cases, namely — 41±6 K, and Mtotai — 65±10 Mq. 
O ur result in dust r nass i s in good agreement with that 
of iZinchenko et all (|2009D . which suggests a total mass 
of 74 M0 derived from dust continuum for W3-SE and a 
virial mass of 90 M©. 

On the other hand, we can estimate the extinction 
from the cold dust component at mid-IR wavelengths 
according to recently studies of the mid-IR extinction 
law. Based on the high-resolution observations with the 
SMA at A 1.1 mm, most of the cold dust (over ~85%) 
is concentrated in the two compact mm/submm sources 
(SMA-1 and -2), the two possible protostellar cores. As- 
suming an even distribution of the dust over the sizes 
of the two cores (see Column 5 of Table 2), the average 
column density of the H2 in the cores would be 1.5x10^"* 
cm-2. On the basis of the equati on [N(HI)-f- 2 N(H2 )]/AK 
= 1.67x10^2 cm-2 mag-i from lLada et all (|2009D . such 
a column density is equivalent to an extinction Ak — 
~200 mag at K band. Assuming a value of 0.5 for both 
AiRAC I -^Ks and Ai/^l Aks b ased on recent studie s on 
the mid-IR extinction law (e.g. lFlahertv et aLll2007f ). the 
corresponding extinctions at the IRAC bands and the 
A 24 /im would be ^100 mag. Obviously the light from 
the sources within the two cold dust cores are severely 
obscured at mid-IR wavelengths. The residual of the 
cold dust emission (~15%, or 0.2 Jy at A 1.1 mm) ap- 
pears to be distributed outside the two compact cores, 
corresponding to a mass of ^^10 M©. On the assump- 
tion of an even distribution of this extended cold dust 
component over a region of 0.15 pcx0.15 pc (15"xl5"), 
the column density N(Il2) is about 2.5x10^^ cm~^, cor- 
responding to an extinction of ^^3 mag at K band, or 
^^1.5 mag at the IRAC bands and A 24 /im. Consider- 
ing the observed filamentary structure of the cold dust in 
the SMA and CARA4A observation, the extinction of the 
extended cold dust distribution to the mid-IR emission 
from W3-SE varies around 1 mag. This extinction may 
have some effects on obscuring the emission from the hot 
dust component depending on the relative distribution 
of the extended cold dust with respected to the hot dust 
component. 

Nevertheless, our fitting to the observed SED of W3-SE 
suggested the existence of a dominant cold dust compo- 
nent with a dust temperature of 41±6 K and a total mass 
of 65±10 M©. A secondary dust component which has 
a little mass but high dust temperature is necessary to 
account for the observed emission at IRAC wavelengths, 
which is likely heated by the nearby stellar objects (e.g. 
IRAC 56). The suggested heating process appears to be 
also consistent with the structure of the mid-IR emission 
characterized as several discrete compact components as 
discussed in the previous section (also see Figure 2). 

3.2. Molecular Lines 
3.2.1. N on- detection of SiO (2-1) Line 

On the basis of CARMA observations of SiO(2-l) line 
at 86.243 GHz, we did not detect SiO emission at a 3cr 
level (^0.3 Jy beam~^) with a channel width of 0.42 
km s^^. The SiO(2-l) line is not a good tracer for the 



outflows in W3-SE. 

3.2.2. Distribution of the HCO+ (1-0) Gas 

Figure 4 shows the channel maps of HCO+(1-0) from 
W3-SE with the CARMA observations, revealing com- 
plex structures in kinematics and spatial distribution. 
Strong line emissions are present in the velocity range 
from -43.9 to -33.2 km s~^ (above 3 tr), with a gap 
between -38.6 and -37.3 km s~^. The emission blue- 
shifted with respect to the gap is stronger and more ex- 
tended than the red-shifted gas. Figure 5(a) shows the 
IICO+(1-0) line intensity image integrated over all the 
line channels, delineating an arc-like morphology in the 
distribution of molecular gas. Hereafter, we denote the 
main component with the maximum line intensity as A, 
the northern and SE lobes as B and C, respectively. The 
main component A is located at the position (-5.4", - 
1.5") offset from the reference center, which coincides 
with the continuum peak position at A 3.4 mm and the 
centroid position of the SMA double source. The spec- 
tra derived from the three positions (A, B, and C) are 
shown in Figure 5(b). The spectrum at A clearly shows 
a double-peak line profile; one is at -40.0 km s~^ (or the 
"blue peak" hereafter) and the other at -36.0 km s~^ (or 
"red peak" hereafter). A spectral dip to ~0 Jy between 
the the blue and red peaks is shown, which is observed 
as an empty gap between -37.3 and -38.6 km s^^ in the 
channel maps (Figure 4). In the velocity range from -34.9 
to -36.9 km s~^, the red-shifted emission peak appears to 
coincide near the position A with the blue-shifted emis- 
sion peak in the velocity range from -39.4 to -40.6 km 
s~^, suggesting that the red- and blue-shifted gas com- 
ponents originate from a common core. At the position 
A, the line peak intensities are 4.3 and 2.1 Jy beam^^ 
for the red- and blue-shifted gas in W3-SE, respectively. 

The position B is located at (-5.1", 9.5"), N of the main 
core A, appeared to be a sub-core. At the position B the 
line emission is much weaker than that at the position A, 
but the line intensity ratio of the blue and red peak is ~2, 
similar to that observed from the position A. Towards 
the position C (-1.1", -6.6"), the spectrum is dominated 
by the blue-shifted line component, while the red-shifted 
line emission appears to be very weak. 

Besides the three relatively compact structures, there 
appears to be an extended halo, ~40" in W3-SE. The 
total HCO+(1-0) flux by integrating the line emission 
above 3ct from the entire W3-SE region is 124 Jy km s~^, 
which underestimates the true flux because of missing 
short spacing. The shortest projected baseline of the 
CARMA data is 2.7 kA, suggesting missing information 
on the source structure with a size larger than ^^76". We 
combined the CARMA data with the IRAM 30-m data 
to assess the effects due to missing short spacing (§2.2). 

Figure 6(a) shows the image constructed from com- 
bined CARMA and IRAM 30-m data with the same 
beam size (6.3" x 5.7", P. A = -81°) as that of the 
CARMA image(Figure 5(a)). The combined image 
shows better the emission from the extended halo but 
trades off the detailed structures. The total line flux 
of 260 Jy km s~^ derived from the combined image is 
about twice as large as that derived from the CARMA 
data alone. Towards the peak position A, the spectrum 
shows that the intensity peak from the combined data is 
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~25% higher than that of the CARMA data while the 
intensity at the spectral dip increases to 0.5 Jy (5 a) in 
the combined image from '^0 Jy of the CARMA image. 

3.2.3. Outflows 

The broad wings of the hne profiles in the HCO+(1-0) 
spectra suggest that outflow(s) is active in this molecular 
cloud. We imaged the high-velocity portions of the out- 
flow(s) by integrating the blue spectral wings from -41.4 
to -46.7 km s~^ and the red spectral wings from -29.6 
to -36.5 km s~^. Figure 7(a) shows the blue- and red- 
shifted gaseous components in blue and red contours, re- 
spectively. It appears that the red-shifted outflow lobe is 
located NW of the continuum core while the blue-shifted 
outflow lobe shows two components located SE (C) and 
N (B) of the continuum core, respectively. The com- 
plex outflow structure implies that there might be more 
than one outflow in W3-SE, and/or the direction of the 
outflow is close to the line of sight. Cutting along a di- 
rection from NW to SE (marked as the straight line with 
P.A. = 120° in figure. 7(a)), we constructed a position- 
velocity (P-V) diagram, showing a kinematic structure 
whic h appears to be consistent with a decelerated out- 
flow (|Qin et al.ll2008l ). The P-V diagram shows only the 
high velocity components that originate from the major 
core A. However, the major axis of the large-scale out- 
flow is poorly determined due to the poor coUimation of 
red- and blue-shifted emission. 

The line emission feature B mimics an independent 
molecular core based on the overlapping blue- and red- 
shifted outflows shown in Figure 7(a) although we do 
not have further evidence either supporting or ruling out 
the possibility of it being a low-mass star forming core. 
The feature B can also be produced from possible multi- 
ple outflows originated from the major core A. The line 
emission feature C appears to be dominated by a blue- 
shifted outflow lobe originated from the major core A. To 
clarify the multiplicity of the outflows in W3-SE, higher- 
resolution observations are necessary. 

3.2.4. Possible Infall Signature 

Blu e-shifted line profiles were discovered by iWu et afl 
(|2007[ ) in HCO+(1-0) and othe r optically thick lines , 
which serve as an infall signature (jMardones et al.lll997f ) . 
Figure 5(b) shows the spectra at three positions of the 
components A, B and C. All the spectral profiles to- 
wards the three regions appear to be characterized by 
an asymmetric line profile with a dominated bl ue-shifted 
peak a nd a weaker red-shifted peak. However, iWu et all 
(|2007f ) also noticed the existence of the significant line 
wings in HCO+ lines from their IRAM 30-m observation. 
From our CARMA observations, the observed asymmet- 
ric double-peak profiles in spectra towards B and C ap- 
pear to be due to outflow. We have shown good evi- 
dence for star formation activities occurring in W3-SE 
and kinematics of the molecular gas in the central 50" 
region which is likely dominated by outflows. We made 
an attempt to investigate the possible contribution from 
infall in the line proflle of the spectrum towards A, the 
dust core, with a model including both infall and outflow. 

In order to quantitatively differentiate the contribu- 
tions from different dynamical processes in the resulting 
spectral proflle, we utilize the infall model of lMvers et al.l 



(|1996f ) to flt the spectral line profile. Four parameters: 
peak optical depths (tq), kinetic temperature (Tk), ve- 
locity dispersion (cr), and the characteristic velocity (Mn) 
with respect to the system velocity of (Ksys), are used in 
the model to characterize two layers of infall gas. We 
found that, indeed, it is difficult to flt the observed spec- 
tral proflle using a model with only infall components 
because of the broad blue- and red-shifted wings produc- 
ing sizable residuals in the fltting. 

Thus, we introduced an additional two-layer outflow 
component with four more characteristic parameters, 
namely, brightness temperature Tb.out, optical depth 
(To.out), characteristic velocity (T4ut) with respect to 
the system, and the velocity dispersion (cout) in the 
Gaussian profile. Similar approach has been considered 
in mod eling the infall motions in NGC 1333 IRAS 4 
(|Attard et al.ll2009() . Here, Fi gure 8 shows the result 
of fitting the observed spectral profile with the modified 
infall model. Assuming Tb.out = 50 K and Vout = 2 km 
s^^ the best fitting gives tq = 4.1 ± 0.4, = 32 ± 5 
K, cr = 0.64 ± 0.04 km s^^, Vin = 0.86 ± 0.06 km s'^, 
Vsys — —38.6 ± 0.1 km s^^ for the infall component (the 
dashed profile in Figure 8) and To,out — 0.2 ± 0.02 and 
Cout — 1-7 ± 0.1 km s"^ for the outflow (dotted profile). 
The spectral dip cannot be produced with a simple super- 
imposition of the double-Gaussian outflow profile (inset 
(a) in Figure 8) to the infall proflle (inset (b) in Figure 
8). The spectral dip seen in the resulting line profile is 
due to the absorption of the cold infall gas located in 
front of the outflow. After solving the radiative transfer 
equations, the final resulting proflle of the model fits the 
data reasonably well. 

The result for the outflow component is not unique. 
Other combinations of the three parameters Tb,out) Kjut 
and To, out are also possible to produce a good fit to the 
observed line profile. In addition to the second kine- 
matic component due to outfiow, any other kinematic 
components within the CARMA beam (~6") character- 
ized with two Gaussians separated in velocity, such as 
two sub-cores orbiting around the mass center, can be 
incorporated into the model in a way similar to that of 
the outflow, i.e., the dynamics in the central 6" may be 
equally important in the production of the observed spec- 
tral profile. Furthermore, we also note that the observed 
asymmetric line profile can be produced by an outfiow 
alone if an unequal optical depth is used for red- and 
blue-shifted components in the outflow. Although the 
modeling is suggestive, we need improved observations 
with high-angular resolution to differentiate between the 
dynamical processes. 

The Vsys = -38.6 km s~^ derived from the model fltting 
differs from the mean value (-38.9 km s~^) of the peak 
velocities from the single-peak opticall y thin lines suc h 
as C^^O (1-0) observed with single-dish (jWu et al.ll2007f l. 
Such a mean peak velocity derived from observations of 
optically thin molecular lines is closer to the systematic 
velocity of the cores. The small discrepancy between the 
values of Ksys predicted by our model and derived from 
observations might be due to the assumption of equal 
optical depth in the blue- and red-shifted outflows used 
in the model. Furthermore, the absorption dip in the 
observed HCO+(1-0) spectrum is centered at ~-38.0 km 
s""'^, red-shifted relative to the systematic velocity (-38.9 
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km s^^) by 0.9 km s~^, which is in good agreement with 
the value of Vin derived from the model fitting. 

With the caveats discussed above, we can calculate the 
kinematic infall rate dM/dt = 47rr^mH/^'^nrVin assuming 
the spherical geometry. Given a radius r = 6x 10'^ AU for 
the main core A, a critical density of HCO"'"(1-0) rim — 
1.8x 10^ cm^'^ (lEvansll999[ ) and a mean molecular weight 
fi = 2.35, an infall rate dM/dt = 9.7x10"^ Mq yi'^ is 
inferred. This infall rate appears to be consistent with 
the global gravitational rate a^/G = 6.8x10"^ Mq yr~^ 
(|Shulll977n assuming that the isothermal sound speed a 
equals to the resulting velocity dispersion cr (0.6 km s~^). 
Indeed, the sound speed of ideal H2 gas at T = 40 K is 
y/1.4:kT/mB_2 —0.5 km s^^, consistent with the value of 
the fitted tr. 

Alternatively, rotation along with other dynamical pro- 
cess and the depletion of HCO+ in the range of veloc- 
ity close to the systematical velocity might be also re- 
sponsib le to the observed asymme tric double-peak line 
profile (jPavlvuchenkov et al.l 120081 ) . Especially, the ve- 
locity gradient as indicated in the P-V diagram o f the 
optically thin 01^0(1-0) line (see iWu et alTl2007[ ) sug- 
gests that a rotation may indeed play an important role 
in formation of the observed line profile. In addition, 
although core A has been resolved into two continuum 
cores in the SMA observations, this multiplicity does not 
account for the observed double-peak line profile since 
further molecular-line observations with higher angular 
resolution show that the double-pea k line profile exists 
in both of the cores ()Zhu et al.ll2009f ). 

4. DISCUSSION 

4.1. The Dynamical Age and Mass of the Outflows 

The dynamical age of the outfiow can be inferred from 
the spatial extension and the mean velocity of the out- 
flow. In the case of W3-SE, the maximum relative ve- 
locities (>5 cr) of the blue- and red-shifted outflow lobes 
from the core A are ~6 km (see the P-V diagram 
in Figure 7(b)) w ith respect to t he systematic velocity 
of -38.9 km s'^ (|Wu et al.ir2007f ). Assuming a size - 
20" (~4xl0^ AU) for each lobe, the dynamical age of 
the outflow is ^^3x10'* yr. Since the major axis of the 
outflow is poorly determined and its direction might be 
close to the line of sight, the estimated dynamical age 
corresponds to a low limit. Given the non-detection of 
radio continuum, both the typical outflow velocity and 
the dynamical age indicate that W3-SE probably hosts 
intermediate- and/or low-mass stars. 

We calculated the mass associated with the main out- 
flow in W3-SE using the formula for column density in 
the case of a hnear, rigid rotor m olecule on the assump- 
tion of LTE (IScoville et all [1981 . For HCO+(1-0), the 
total column density at all energy levels can be estimated 
following the equation: 

(1 - e^T^t)-^ J Tydv (3) 

where we adopted the rotational constant B = 44.594 
GHz and the permanent electric dipole moment — 
3.30 debye ifWoods et al..,1975.) . Assuming optically thin 



gas and excitation temperature of Tqx = 40 K, the optical 
depth r is proportio nal to the brightness temperature Tb 
(|Wilson et al.ll2009D : 

^ kTs , 1 1 

hu e'"'/'=^" - 1 e'^"/''^''!! - 1 

The outflow mass can be calculated from Equation (3) 
in the case of optically thin emission, which valid for the 
outflow gas in the high- velocity wings, 

M{ll2) = X^l,^^m{B.2)D^ X 

j iV(HCO+)df] (4) 

where we adopted a fractional abundance Xhco+ — 
7.6x10^9 for W3-SE (Zinchenko et al . 2009) and a dis- 
tance D = 2 kpc. Then, we determined the velocity- 
integrated brightness temperature averaged over the 
emission region ~ 50"x50", J < Tb >n'iv = 1.2 K km 
between -46.7 and -41.8 km s~\ and J < Tb >ndv 
= 1.4 K km s~^ between -36.5 and -31.6 km s~^. 
The corresponding masses are 2.5 and 2.8 Mq for the 
blue- and red-shifted outflow lobes, respectively. Conse- 
quently, with the maximum relative velocities of the out- 
flow lobes (~6 km s^^) and the dynamic age (~3xl0''^ 
yr) , the momentum rates for both the blue and red lobes 
are ~5xl0^^ M0 km s^^ yr^^, corresponding to me- 
chanical luminosities of 0.2 Lq. The inferred outflow 
momentum rates and mechanical luminositie s appear to 
be co nsistent with B-type stars in W3-SE (jArce et al.l 
I2OOI . 

4.2. A Group of Young Stars /Protostars in W3-SE 

Observations at both millimeter and submm wave- 
lengths, together with those at mid-IR suggest that the 
main molecular core A hosts a group of young stars 
and/or protostars. The millimeter and submm observa- 
tions suggest that the cold dust is distributed in an elon- 
gated region from SE to NW (Figure 2) with two dust 
condensations separated by 4". The non- detection of ra- 
dio continuum at A 6 cm (|Tieftrunk et al.iri998) indicates 
that no OB stars have been found in W3-SE yet although 
we cannot rule out the possibility that the two sub-cores 
in the molecular core A (SMA-1 and -2) will continue to 
accrete from the molecular gas reservoir of ^65 M0 and 
form a massive star. The four IRAC images show a small 
group of discrete sources in the dust core. One (I RAC 
56) of them has been identifled bv lRuch et al.l ()2007D as a 
point source, suggesting an association with a YSO ear- 
lier than Class I. This source appears to be one of the 
youngest object in the entire W 3 clus t er as indicated in 
the color-color diagram of Ruch et all ()2007n . The weak 
point-like source IRAC 56a might be associated with a 
protostar embedded in the dust core SMA-1. 

The positions of the emission peaks in the millime- 
ter/submillimeter continuum and lines appear to be dis- 
placed significantly from the brightest mid-IR sources in 
this region. The majority of the mid-IR emission arises 
from the region located east of the molecular core (see 
Figure 2). The displacement likely indicates that the 
hot and cold dust components which are inferred from 
fitting the SED are also spatially separated in the W3- 
SE region. The spatial separation between the hot and 
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cold dust might be explained by sequential star forma- 
tion in W3-SE. Namely, stars formed first in the region 
east to the main molecular core A in W3-SE, where the 
dust/molecular components might have been substan- 
tially dissociated or destroyed by the stars. The activ- 
ity of forming new stars is now taking place in A. The 
activity of star formation appears to propagate from E 
to W in W3-SE. Alternatively, the displacement of the 
emission peak position at mid-IR from that of millime- 
ter/submillimeter could be due to possible strong mid-IR 
extinction in the main molecular core A. The scenario of 
sequential activities in star formation needs to be veri- 
fied by further improved mid-IR observations with cor- 
rections for dust extinction. 

5. SUMMARY AND CONCLUSIONS 

We observed W3-SE with CARMA and SMA at A 3.4 
and 1.1 mm. Using the continuum flux densities deter- 
mined with CARMA, SMA and JCMT together with 
measurements at mid-IR wavelengths with Spitzer and 
MSX, we constructed the SED for the W3-SE core. The 
SED can be fitted with two thermal dust components, 
suggesting that the majority of the gas/dust (65±10 Mq) 
is at a relatively low temperature of 41±6 K correspond- 
ing to radiation peaked at a far-IR wavelength. There is 
an additional hot component (~400 K) with little mass 
accounting for the secondary emission bump at mid-IR 
in the SED, suggesting that a small fraction of dust has 
been heated by newly-born stars in the immediate sur- 
roundings. 

The CARMA observations of the HCO+(1-0) line con- 
firmed the presence of blue-shifted double-pea k line pro- 
file to wards the W3-SE core discovered by iWu et ahl 
(|2007f) with the IRAM-30m. We made an attempt to 
fit the line profile with the infall model of Myer s et ahl 
()19 96^ by incorporating outflow components. The model 
fitting suggests that in addition to infalling gas, a signif- 
icant contribution from outflows is necessary to produce 
the observed line proflles. The kinematics of the molec- 
ular gas in W3-SE observed with CARMA in the central 
40" region are consistent with a bipolar outflow with a 
dynamic age of ~3xl0'' yr. Finally, based on the ob- 
servations with CARMA, SMA, and JCMT along with 
the previous published results, we conclude that W3- 
SE hosts a group of young stars and protostars which 
are likely located in the double dust cores observed with 
SMA. 
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TABLE 1 

Calibrators for the CARMA and SMA Observations 



Date 


A 

(mm) 


Bandpass 


Phase 




Flux Density 


yource 


Flux(Jy) 


yource 


HuxlJy) 


Source 


i^'lux(Jy) 


2008-06-24 


3.4 


3C84 


6.1 


0102-1-584 


2.1 


Uranus 


7.2 


2008-06-27 


3.4 


3C454.3 


27.8 


0102-1-584 


2.5 


MWC349 


1.3 


2008-10-27 


1.1 


3C454.3 


13.2 


0244-1-624 


0.6 


Uranus 


45.6 



TABLE 2 

Measurements of the Continuum at mm/submm 



A 

(mm) 


Telescope 




Gaussian Fitting Results 




(Jy) 


Peak(") 


Ipeak(Jy beam ^) 


Sizes(") 


P.A.(°) 


3.4 


CARMA 


-5.5,-0.9 


0.024±0.002 


7.9x5.9 


-51 


0.04±0.01 


1.1 


yMA 










1.4±0.4 




(SMA-1) 


-6.2, 0.2 


0.21±0.01 


4.4x3.0 


21 


0.72±0.01 




(SMA-2) 


-4.2,-2.8 


0.17±0.01 


3.5x1.5 


76 


0.35±0.01 


0.85 


JCMT 


-8.2, 1.2 


3.3±0.1 


25x18 


-48 


8±1 


0.45 


JCMT 


-8.8, 0.5 


17.1±0.7 


18x11 


-41 


75±8 



Note. — Columns 1 and 2 give the wavelengths and telescopes which are used in the observations. Columns 3, 4, 5 and 6 are the results 
from Gaussian fitting to the sources, corresponding to peak positions, peak intensities, core sizes and position angles. Column 5 gives the 
total flux densities integrating the emission >5 cr over the continuum core region of 15"xl5" (except for SMA-1 and -2). The uncertainties 
are attributed mainly to the missing short spacing data for the interferometer observations and the uncertainty in determination of the 
boundary for the continuum core in the JCMT images. For SMA-1 and SMA-2, the total flux densities are determined from multiple 
Gaussian fitting to the 1.1-mm dust emission. 



TABLE 3 

The Continuum Flux Densities at mid-IR 



A 


Instrument 


S,. 


(/im) 




(Jy) 


24 


MIPS 


8.8±1.3 


21.3 


MSX 


9.2±3.7 


14.7 


Myx 


0.83±0.33 


12.1 


MSX 


0.94±0.38 


8.3 


MSX 


0.95±0.38 


8.0 


IRAC 


0.82±0.16 


5.8 


IRAC 


0.48±0.10 


4.5 


IRAC 


0.18±0.04 


3.6 


IRAC 


0.054±0.011 



Note. — The observing wavelengths, the relevant instrument and the continuum flux densities are given in columns 1, 2 and 3, 
respectively. 
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Fig. 1. — The dust emission of W3-SE. (a) The image of the continuum emission at A 3.4 mm made from the CARMA observations by 
combining two 500-MHz bands. The contours are Scrx (-1, 1, 2, 3, ...) and 1 tr = 0.0015 Jy beam~^. The FWHM beam size given at 
bottom-right is 6.2"x5.6", P.A.=-81°. (b) The image of continuum emission at 1.1-mm based on the SMA observations in the compact 
array with two 2-GHz side bands. The contours are 5(7X (-1, 1, 2, 3, ...) and 1 tr = 0.0015 Jy beam"!. The FWHM beam size is 2.9" X 2.7", 
P. A. = —64°. The open circle in both (a) and (b) marks the peak position of the 3.4-mm emission. The images of continuum emission at 
0.85 mm (c) and 0.45 mm (d) were made with the JCMT observations with the FWHM beam sizes of 14.5" and 7.5", respectively. The 
crosses in (b), (c) and (d) mark the peak positions of the double sources (SMA-1 and -2). The units of the wedge in (a) and (b) are Jy 
beam~^; mjy beam~^ for (c) and (d). 
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Fig. 2. — Spitzer images at (a) MIPS 24 /im and (b) IRAC 3.6 fim. In the two images the wedge scales the continuum emission in units of 
10^ Jy sr~^, which is equivalent to 2.35 X 10~^ Jy arcsec"^ . The crosses in the two images indicate the peak positions of t he double so urces. 
The filled circle in (a) and (b) marks position of the mid-IR point source (IRAC 56) which has been identified bv lRuch et al.l ((2003). The 
open circles in (b) indicate two additional sources (IRAC 56a and 56b) in the central core region of W3-SE. 



12 



Zhu et al. 




Fig. 3. — The SED of the continuum emission from the W3-SE core determined from the observations with the CARMA, SMA, JCMT, 
MSX, MIPS and IRAC. The legend (top-right) shows the symbols corresponding to each of the telescopes. The thick solid curve outlines 
the resulting SED fitting from two thermal dust components (indicated by the thin curves) determined from our model fitting, (a) The 
best (LSQ) fit to the flux density measurements of the SED assuming the optically thin dust emission from both cold and hot components, 
(b) A good fit to the SED considering the effect of self-absorption of the cold component with /3 = 2 and assuming that the cold component 
does not overlap the hot dust emission region. The turnover frequency J't=i of the cold component is shifted by Ai/ = 80 GHz to higher 
frequency as compared to that derived from the optically-thin approach. 
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Fig. 4. — Channel maps of the HCO+(1-0) line observed with the CARMA. The channel separation is 0.42 km s ^. The contours 
are 4(tx(-1, 1, 2, 3, ...) and 1 a = 0.1 Jy beam"'^. The LSR velocity is labeled at top-right on each channel map. The FWHM beam 
(bottom-right) is 6.3"x5.7", P.A.=-79°. 
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Fig. 5. — (a) The integrated intensity map of HCO+(1-0) from the CARMA observations. The eontours are 3(tx(-1, 1, 2, 3, ...) and 1 
a = 0.6 Jy beam^^ km s~^. The FWHM beam size is 6.3"x5.7", P.A. = -81°, shown at bottom-right. The eross marks the A 3.4-mm 
continuum peak, (b) The spectra towards positions A, B and C marked in Figure 5(a). 
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Fig. 6. — (a) The integrated intensity image of the HCO^(l-O) line emission constructed by combining the CARMA data with the IRAM 
30-m data. The beam size of 6.3" x5. 7", P.A. =-81° is the same as the one of CARMA image (Figure 5(a)). The boundary eflfect has been 
discussed in the text, (b) The top spectrum is the integrated HCO+(1-0) line emission within the 25" X 25" square as delineated. The 
spectrum toward the peak position A is shown at the bottom. 
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Fig. 7. — The integrated intensity maps for the blue and red wings based on CARMA data. The contours are 5<jx{1, 2, 3, ...) and I a = 
0.12 Jy beam-i km s'^. The FWHM beam size is 6.3"x5.7", P.A.=-81°. The straight Hne denotes the direction of the cut for the P-V 
diagram in (b). (b) The P-V diagram along the cut indicated in Figure 7(a). The vertical axis is the offset relative to the position A, in 
units of arcsec, positive to SE. The contours are 5crx(l, 2, 3, ...) and 1 a = 0.08 Jy beam^^. 
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Fig. 8. — The observed line profile (crosses connected with a black curve) is fitted with a two-layer infall model including blue- and 
red-shifted outflow components. The solid curve (red) shows the resulting profile from the infall (dashed curve) and outflow (dotted curve) 
components. The insets at the top-right show the two spectra produced separately with infall alone (upper) and outflow alone (lower) 
using the same parameters derived from the overall spectra fltting. 



